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The observation of extensive fractional quantum Hall states in graphene brings out the possibility 
of more accurate quantitative comparisons between theory and experiment than previously possible, 
because of the negligibility of finite width corrections. We obtain accurate phase diagram for dif¬ 
ferently spin-polarized fractional quantum Hall states, and also estimate the effect of Landau level 
mixing using the modified interaction pseudopotentials given in the literature. We find that the 
observed phase diagram is in good quantitative agreement with theory that neglects Landau level 
mixing, but the agreement becomes significantly worse when Landau level mixing is incorporated 
assuming that the corrections to the energies are linear in the Landau level mixing parameter A. 

This implies that a first order perturbation theory in A is inadequate for the current experimental 
systems, for which A is typically on the order of or greater than one. We also test the accuracy of 
the composite-fermion theory and find that all lowest Landau level projection methods used in the 
literature are very accurate for the states of the form nf{2n + 1) but for the states at nf{2n — 1) 
the results are more sensitive to the projection method. An earlier prediction of an absence of spin 
transitions for the n/(4n+ 1) states is confirmed by more rigorous calculations, and new predictions 
are made regarding spin physics for the nf{4n — 1) states. 

PACS numbers: 73.43.Cd, 71.10.Pm 


I. INTRODUCTION 

A longstanding issue in the field of the fractional quan¬ 
tum Hall effect (FQHE]P has been that the quantitative 
agreement between theory and experiment is less precise 
than what one would expect. It is possible to obtain 
very accurate numbers for many quantities of interest 
from the composite fermion (CF) theor}E®l. Detailed 
comparisons have been carried out for activations gap^, 
collective mode dispersion^, and spin-polarization phase 
transitions^. In all cases, the measured numbers are gen¬ 
erally consistent with those predicted by theory, but the 
agreement is worse than that suggested by the accuracy 
of the theory as determined from comparisons with exact 
diagonalization resultS^^^K^. It is believed that the devi¬ 
ation arises from corrections due to effects extraneous to 
the FQHE physics, such as finite quantum well width, 
Landau level (LL) mixing and disorder, which provide 
sizable corrections that are hard to deal with in a quan¬ 
titative manner. 

The observation of FQHE in graphend^^^^ provides a 
unique opportunity in this context, because finite width 
corrections are essentially absent in graphene. Further¬ 
more, as noted by Peterson and Nayakl^, unlike in GaAs 
quantum wells, in the n = 0 LL of graphene, LL mix¬ 
ing does not produce any effective three body interac¬ 
tion (which incorporates the breaking of particle-hole 
symmetry), but only corrections to the pairwise inter¬ 
action. One might therefore expect smaller corrections 
due to LL mixing in graphene than in GaAs. The FQHE 


in graphene may thus provide an opportunity for bet¬ 
ter quantitative comparisons between theory and exper¬ 
iments, and a better appreciation of our understanding 
of the role of LL mixing. 

The spin polarization transitions provide some of the 
most precise tests of the quantitative accuracy of the 
FQHE theory, for several reasons. First, it is a ther¬ 
modynamic measurement (as opposed to the excitation 
energies of charged or neutral modes), and therefore is 
likely to be less susceptible to the presence of disorder. 
Second, the critical Zeeman energy where a transition be¬ 
tween two differently spin polarized FQHE states occurs 
provides a direct measure of the rather small Coulomb en¬ 
ergy differences between the two states, and thus enables 
a detailed and exceedingly sensitive test of our quanti¬ 
tative understanding of the FQHE. Finally, there is an 
extensive amount of experimental phenomenology asso¬ 
ciated with the spin physics. Phase transitions as a func¬ 
tion of the Zeeman energy have been measured in various 
semiconductor based two-dimensional system^^^M^, and 
recently also in graphene by Feldman et al]^. 

Hoping to obtain a better comparison between the¬ 
ory and experiment, we have determined the spin- 
polarization phase diagram as accurately as possible for 
an ideal two-dimensional system in the n = 0 LL. (These 
results apply to graphene as well as to narrow GaAs 
quantum wells, because in the limit of zero thickness, 
the n = 0 LL wave functions are identical for the two.) 
We perform large scale exact diagonalization studies of 
various FQHE states for this purpose. This also allows 
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us to determi ne the quantitative accuracy of the CF 
theor3^ ^F4 | 35 [ estimate the effect of LL mixing 

by modifying the pseudopotentials according to RefM 
Our conclusions, briefly, are as follows. 

From exact diagonalization calculations, we have ob¬ 
tained essentially the exact theoretical phase diagram of 
the spin polarization of FQHE states without allowing for 
LL mixing (see blue crosses in Fig.[^. We find that it is in 
excellent agreement with the experimental phase diagram 
in graphene (see dots in Fig. [^. The agreement between 
theory and experiment becomes significantly worse when 
LL mixing is included to linear order in the LL mixing 
parameter A using the pseudopotentials of Ref. |T7| This 
suggests that the current experimental systems are out¬ 
side the linear regime, which is not surprising given that 
A is typically of order 1 and sometimes much larger (e.g. 
2.2 for suspended graphene). This is relevant to the issue 
of whether the Pfafiian or the anti-Pfafiian stat^^^M^ at 
z/ = 1/2 in the second LL is selected by LL mixing, which 
has received much attention recent 

The allowed spin polarizations for various FQHE states 
at n/(2n ± 1) and their energy ordering are correctly 
predicted by the CE theory. The mean-field model in 
which composite fermions a re trea ted as free particles 
at an effective magnetic fielcP^^^^^ remains satisfactory, 
and we obtain a precise estimation for the CE mass. We 
also carry out quantitative tests of the CE theory, which 
are extremely precise because the critical Zeeman ener¬ 
gies depend on very small energy differences between dif¬ 
ferent states, and even a slight error in the energy can 
lead to large corrections in the critical Zeeman energies. 
Eor the fractions n = njifln + 1), calculations based on 
Jain’s wave functional predict the critical Zeeman ener¬ 
gies with 15% accuracy, which can be further improved 
by the method of CE diagonalization that incorporates 
A level (AL) mixing, where a AL refers to a CE Lan¬ 
dau level in an effective magnetic field. Eor the states at 
u = nf (2n — 1), the results depend sensitively on the low¬ 
est Landau level (LLL) projection method. In part icular, 
the treatment with Jain-Kamilla (JK) projectiorPE^El 
underestimates the critical Zeeman energies by a factor 
of 2-3 (i.e., relatively overestimates the energies of the 
non-fully spin polarized states). This is not an intrinsic 
deficiency of the CE theory but rather a technical issue, 
as can be seen from the fact that the “hard-core pro¬ 
jection” introduced previously in Ref. [8] produces very 
accurate states for n = njifln — 1). (This projection is 
presently not amenable to calculations for large systems.) 

Eor completeness, we have also considered the states of 
composite fermions with four vortices attached. A previ¬ 
ous predictioiP^ of the absence of spin transitions for the 
n/(4n + 1) states is confirmed by more rigorous calcula¬ 
tions. Here the states remain spontaneously spin polar¬ 
ized even in the absence of a Zeeman energy. Spin tran¬ 
sitions are, however, possible for the ^CE EQHE state 
at z/ = n/(4n — 1), and we estimate the critical Zeeman 
energies for the prominent transition. 


II. FQHE IN GRAPHENE 

The physics of graphene differs from that of GaAs in 
two important aspects: the dispersion of electrons is lin¬ 
ear and there are four Dirac cones. The linear dispersion 
of Dirac fermions leads to LLs which have a cyclotron en¬ 
ergy of ^gIi{n)^J2\n\hvF /where Ui? is the Fermi veloc¬ 
ity and ^ is the magnetic length and n is any integei®^. 
(This is to be contrasted from GaAs or other conven¬ 
tional semiconductors wherein electrons have a parabolic 
dispersion and LLs which have a cyclotron energy given 
by hujc{n + 1/2), where n is a non-negative integer and 
ujc is the cyclotron energy. ) Also, the LL wave func¬ 
tions are two component wave functions, corresponding 
to two sublattices of graphene. In general this leads to 
different interaction pseudopotential^^^^^ than in sys¬ 
tems with parabolic dispersion. However, it turns out 
that the electron wave functions as well as the interac¬ 
tion pseudopotentials in the n = 0 LL of graphene are 
identical to those in the LLL of GaAs. 

The second difference is that each LL of graphene has 
four bands, which arise from the valley and spin degrees 
of freedom, while GaAs has two bands from the spin de¬ 
gree of freedom. Within each band, the physics of FQHE 
in the n = 0 LL is identical in the two systems, apart from 
corrections due to finite width and LL mixing. Much 
work has been done toward understanding the origin of 
the lifting of various degeneracie^^^H^, but we will as¬ 
sume below that all bands are well separated. This as¬ 
sumption is justified for the experiments with which we 
compare our results. 

We refer the reader interested in further information to 
an extensive liter ature on the physics of Dirac fermions 
in a magnetic fielcP^^^. 


HI. PHASE DIAGRAM OF SPINFUL CF 
STATES 

A good qualitative and semiquant it at ive theoretical 
understanding of these transitions has been obtained in 
terms of integer or fract ional qu antum Hall effect of 
spinful composite fermionJ- ^F4 | 47 | 57tf 65]^ which successfully 
predicts the allowed spin polarizations at all of these 
filling factors and also provides an estimate of the the 
critical Zeeman energy where transitions between them 
occur. While these quantitative estimates are a good ze¬ 
roth order approximation, their accuracy has not been 
carefully evaluated in the past. 

The FQHE state at z/ = n/ (2pn ± 1) maps into integer 
quantum Hall effect (IQHE) state of composite fermions 
with n filled A levels, where a composite fermion is de¬ 
fined as a bound state of an electron and 2p vortices. 
Eor spinful composite fermions, the CE filling is writ¬ 
ten as n = where and are the num¬ 

ber of filled spin-up and spin-down ALs. The different 
states will be denoted as (n>^,n^), and we will use the 
convention without loss of generality. One 
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FIG. 1: (Color online) Spin phase diagram showing the critical Zeeman energies for transitions of the fractional quantum Hall 
states at the filling factor v — nl{2n zb 1). We define k = EzjEc, where Ez is the Zeeman splitting and Ec = /el is the 
Coulomb energy scale. The crosses show the results from exact diagonalization (assuming zero thickness), the dots from JK 
wave functions, and squares are fits from the free-CF model (see text). The blue, red and green colors indicate transition from 
(n, 0) ^ (n — 1,1), (n — 1,1) ^ (n — 2, 2) and (n — 2,2) ^ (n — 3, 3) respectively. We note that the exact diagonalization results 
for the spin transitions at v — and z/ = 4/9 are obtained using results extrapolated from only two finite systems, and may 
therefore be less accurate. 
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can list all possible states and their spin polarization 
7 = {n^ — n^)/(n^ + n^). For example, for 2/5 and 2/3, 
which both map into n = 2 of composite fermions, we 
have a fully spin polarized state (2,0) (with 7 = 1 ) and 
a spin singlet state (1,1) (with 7 = 0 ). To take another 
example, 6/13 and 6/11 map into n = 6, where we have 
four possible states (6,0), (5,1), (4,2) and (3,3), with 
7 = 1, 2/3, 1/3, and 0, respectively. One expects one 
transition at 2/5 and 2/3 and three at 6/13 and 6/11. 
The possible states and spin polarizations of other frac¬ 
tions can be similarly enumerated. 

The CF theory also identifies the flux values where 
these states occur in the spherical geometr}!^. These are 
the flux values at which our calculations below are carried 
out. All our calculations are performed in the spherical 
geometry, where the ground states have total orbital an¬ 
gular momentum L = 0 and a total spin S that can be 
ascertained by the CF theory. We will assume the ideal 
limit of zero thickness and neglect disorder. In quoting 
the energies below, we include the electron-background 
and background-background interaction. The density for 
a finite system (in the spherical geometry) depends on the 
number of particles N and is slightly different from its 
thermodynamic limit. To eliminate this effect we use the 
density corrected energ 3 J^F^^ = {2Qiy/Ny^‘^EN for ex¬ 
trapolation to the thermodynamic limit N~^ 0, where 

the integer 2Q is the magnetic flux (in units of 0o = he/e) 
to which the electrons are subjected. All energies quoted 


below are the thermodynamic limits of the per particle 
density corrected energies lim7v^c»^Ar/^- 

To avoid clutter, we give only give the spin-polarization 
phase diagram in the main text. All of the individual 
numbers as well as extrapolations are given in various 
tables and figures in Appendix [A} 


A. Exact diagonalization 


We first obtain the extrapolated values of energies of 
the variously polarized states at fractions of the form 
n/(2n ± 1). These include the largest systems for which 
exact diagonalization can currently be performed (see 
Appendix for Hilbert space dimensions). For filling 
factors 4/7 and 4/9, the extrapolated values are obtained 
with only two points and thus must be treated with cau¬ 
tion, but we have chosen to include them because linear 
extrapolation in 1/A has been found to be quite accurate 
for other systems for which several points are available. 
Once we have the energies, we obtain the critical Zee- 
man energies by setting the energy difference of the two 
successive states (n>^,n^) and + 1) to zero: 




( 1 ) 


where we have used fo denote the thermodynamic 

limit of the per particle density-corrected, background- 
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subtracted Coulomb interaction energy of the state 
{n^,nx). This gives 

^ /{etj ~ “ ^(nt-i,n^+i)] (2) 

where e is the dielectric constant of the host material 
and I = ^JhcleB is the magnetic length. The resulting 
critical energies are shown by the blue crosses in Fig. 


B. Free CF model 


We ask to what extent the results may be interpreted 
in a model that treats composite fermions as free particles 
with an effective masJ^, which has bee n used routinely 
to analyze the experimental dat al^^ * ^^ -^. The interaction 
energy between electrons is modeled in terms of the CF 
cyclotron energy, defined as 


e5* 

hujl = h — = h 


eB 


{2pn ± l)m*c 2pn ± 1 


b (3) 


where m* is the CF mass. (The subscript is to remind 
us that this mass is the “polarization mass” of compos¬ 
ite fermions,!^ which is the relevant mass for the spin- 
polarization phase transitions. It is to be distinguished 
from the mass defined from the activation gapP^.) The 
CF mass is often quoted in units of the electron mass in 
vacuum, ruei 

rrie a e^jel 

where Oc = ^Bjm^c is the cyclotron frequency of elec¬ 
tron in vacuum. The CF mass behaves as rrip ^ ^/B and 
for the parameters of GaAs, we have 


m;/me = (0.026/a) (5) 

An immediate qualitative prediction of the free-CF 
model is that the interaction energies increase with the 
degree of spin polarization. This has been found to be the 
case for all states of the form n/(2n± 1). [As seen below, 
this is not the case for n/(4n + 1).] At a more quanti¬ 
tative level, the free-CF model predicts that the critical 
Zeeman energy for the transition between (np,np) and 
(np — 1, np -h 1) is given by 


-[% - % -1]?)^ 


* 

c 


which gives 


—Ez = 0 

n 


K = a 


np - np - 1 

2pn ± 1 


(6) 

( 7 ) 


We have found that the best fit for the critical Zeeman 
energies calculated here is provided by <an/( 2 n+i) = 0.056 
and (an/( 2 n-i) = 0.050. A single value of a gives a slightly 
less satisfactory fit, implying a weak filling factor depen¬ 
dence for the CF mass. Nonetheless, the free-CF model 
works reasonably well. We note that the CF polarization 
mass m* for the reverse flux attached states is about 10% 
higher than that for parallel flux attached states. 


C. Microscopic theory 

The wave function for the (np, np) state at n/ (2pn± 1) 
state is given by 

= VLLL^nJ^^ = (8) 

and 

= ’PLLLi^nVJ^^ = (9) 

where 

j= n (10) 

^<j<k<N 

is the Jastrow factor, Zi is the coordinate of the elec¬ 
tron, (^n^) is the Slater determinant wave function 
for Up (np) filled LLs of electrons, and T^lll denotes LLL 
projection. In this section we shall restrict ourselves to 
the case of p = 1 and consider higher values of p in the 
subsequent section. Three schemes have been employed 
for LLL projection, which result in slightly different LLL 
projected wave functions, (i) “Direct projection” will re¬ 
fer to the method considered in Refs. [7| and [H wherein 
the product wave function is expanded into Slater de¬ 
terminant basis and only the part strictly in the LLL is 
retained. This method can be implemented for relatively 
small systems (ten particles or fewer), (ii) In the “hard¬ 
core projection” of Ref. [H one writes the wave function 
as 


^h^-core ^ j p ^ J (11) 

2n + l 


«'h_!y-co-'e = JPLLL[i>„]*J (12) 

2n—l 

As the name implies, this method explicitly builds corre¬ 
lations such that the wave function vanishes even when 
particles of opposite spin coincide. This method also re¬ 
lies on expansion into Slater determinant basis and can 
be implemented only for small systems, (hi) The most 
widely used projection method is the so-called JK pro¬ 
jection, which has the advantage that it does not require 
expansion into Slater determinant basis and thus can be 
evaluated for very large systems. This method has been 
used extensively to make quantitative predictions for var¬ 
ious quantities. The details of the J K proje ction method 
have been outlined in the literatur^^IMl ^^d will not 
be repeated here. We will refer to the resulting wave 
functions as JK wave functions, to distinguis h them from 
wave functions obtained by other project ion^^HMlI, 

The results from the JK projection are also shown 
in Fig. [2 For the “parallel flux-attached” states at 
n/ (2n + l), the JK wave functions underestimate the crit¬ 
ical Zeeman energies (also given previously in Ref. 
by ^ 15%. As a result, the a ^ 0.056 is lower, and the 
CF effective mass is higher, by ^15 % compared to the 
values from exact diagonalization. For the “reverse-flux 

















5 


ly 

system 

% error in the Coulomb interaction energy 


N 

2Q 

state 

JK 

Direct 

hard-core 

2/3 

8 

11 

(1,1) 

0.75 

- 

0.04 (Ref. [S]) 

2/3 

8 

12 

(2,0) 

0.13 

0.02 (Ref. [H]) 

- 

2/5 

8 

16 

(2,0) 

0.01 

0.01 (Ref. [SZ]) 

- 

3/5 

8 

13 

(2,1) 

0.29 

- 

0.05 (Ref. 0) 

3/5 

9 

16 

(3,0) 

0.03 

0.01 (Ref. 0) 

- 


TABLE I: Percent error in the Coulomb interaction energy 
of the CF wave function obtained from different projection 
methods for several hnite systems for which the exact results 
are known. The last two columns give results from Direct 
projection for fully spin polarized states and from hard-core 
projection for the non-fully spin polarized state. (For the 
fully spin polarized states, Direct and hard-core projections 
give essentially the same energieP^.) The systems correspond 
to N particles on the surface of a sphere subjected to a total 
flux of 2Qhcle. The filling factors and the state are shown. 


attached” states at nl{2n — 1), the JK wave functions 
underestimate the critical Zeeman energies by a factor of 
two to three. These results bring out the limitations of 
the JK projection method for the reverse-flux attached 
states. 

From our current study as well as the previous 
result^^ (reproduced in Tablewe find that: the hard¬ 
core projection produces very accurate energies for fully 
as well as non-fully spin polarized states at n/(2n+l) and 
n/(2n — 1); the JK and the Direct projections produce 
accurate results for fully polarized states at n/(2n + 1) 
and n/(2n — 1) and also for non-fully polarized states at 
n/(2n + 1); the JK and Direct projections are somewhat 
less accurate for non-fully polarized states at n/{2n — 1). 
It is easy to see why the JK / Direct projection under¬ 
estimates the critical Zeeman energies: its deficiency is 
that it does not eliminate configurations in which spin-up 
and spin-down particles are coincident, and thus overes¬ 
timates the energies of non-fully spin polarized states by 
a larger amount than for the fully spin polarized states, 
thus resulting in an underestimation of the critical Zee- 
man energies. It is unclear why the JK and Direct pro¬ 
jections work better for the non-fully polarized states at 
n/ (2n + 1) than those at n/{2n — 1). 

It is in principle possible to improve the accuracy of the 
results within the JK projection scheme by the method 
of composite fermion diagonalization (CFD]I^, in which 
one can obtain more accurate energies by allowing some 
A level mixing. (AL mixing is to be distinguished from 
LL mixing.) We allow A level mixing by including CF 
excitons in the basis. A CF exciton is defined as a pair 
of CF particle and CF hole, where a CF particle is a 
CF in the lowest unoccupied A level and a CF hole is a 
missing CF from the top most occupied A level. Note 
that a single CF exciton does not change the ground 
state energy. This is easy to see since a CF hole car¬ 
ries an orbital angular momentum one smaller than the 
CF particle with which it forms the exciton. Therefore 
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FIG. 2: Schematic representation of the CF basis functions 
used for CF diagonalization study of the 2/(4p + 1) spin- 
singlet FQHE state. Panel (a) shows the “unperturbed” state. 
At the first order approximation, CF diagonalization allows 
hybridization of the state in (a) with the states in (b), (c) 
and (d) to obtain a new ground state with lower energy than 
the unperturbed state in (a). Successive mixing with higher 
and higher excitations produces better approximations to the 
ground state. 


the smallest angular momentum that a CF exciton can 
have is 1, thus precluding its admixture with the ground 
state. Therefore, we need a minimum of two CF exci¬ 
tons to improve the ground state. For example in Fig. 
I^we show the excitations we considered to improve the 
u = 2/(4p+ 1) spin-singlet state. The Hilbert space 
grows very quickly with the number of excitons included 
in the basis for CFD, so we restrict ourselves to at most 
two excitons. Among the wave functions shown above, 
the fully polarized ones are extremely accurate, so this 
procedure of including two CF excitons in the basis of 
CFD only marginally improves the ground state energy 
of the fully polarized state. However for the unpolar¬ 
ized states the improvement is substantial, evidenced by 
the fact that for the spin-singlet states the energies im¬ 
prove by around 10%. The method of CFD produces a 
critical Zeeman energy that is within ^3% of the exact 
value. We have carried out similar first order CFD cal¬ 
culations for many of the states considered in this work. 
These results are shown in the Appendix and are la¬ 
beled as “CFD”. We have found CF diagonalization to 
be impractical for the reverse-flux attached n/{2pn — 1) 
states for technical reasons. (For states involving reverse 
flux attachment, the projected state is obtained as an 
alternating sum of elementary symmetric polynomials of 
high degree; see Ref. |62l To avoid the loss of significant 
digits we have to use software emulated multiple preci¬ 
sion floating point numbers. The number of Monte Carlo 
steps to obtain the overlap and interaction matrices for 
CFD with reasonable accuracy is beyond our reach.) 


IV. COMPARISON WITH EXPERIMENTS IN 
GRAPHENE 

Fig-i shows a comparison between the critical values 
of K obtained from exact diagonalization (blue crosses) 
with the measured ones. The results for graphene are 
taken from Feldman et alM The value of k depends 
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on the product eg where e is the dielectric function and 
g is the Lande g factor, for which we have taken the 
values e = 3.0 and g = 2.0 for graphene in the figure. 
An excellent agreement between the measured and the 
theoretical values is evident. We stress, however, that 
we have not yet included effects of LL mixing, which is 
discussed below. 

For completeness, Fig.j^also shows critical Zeeman en¬ 
ergies measured in other systems, taken from Padman- 
abhan et al. (AlAs quantum well]P^ and Du et. aP^ 
(GaAs-AlGaAs heterojunction). In all cases, the critical 
values of n as well as their filling factor dependences are 
roughly consistent with theory. The difference between 
the critical values of n in graphene and in hetero junction 
samples is somewhat surprising, because the hetero junc¬ 
tion samples also correspond to a very small thicknesses. 
We believe that part of the difference might result from 
the fact that Feldman et al varied k by changing the 
density, whereas Du et al^ accomplished that by tilting 
the magnetic field. We believe that high parallel fields in 
the latter experiments may cause additional corrections 
which have not been considered. 


A. corrections due to LL mixing in graphene 

We next come to the effect of LL mixing. The strength 
of LL mixing is measured by a parameter A, which is 
defined as the ratio of the Goulomb to the cyclotron 
energ}!^. For massive electrons (e.g., in GaAs) A = 
(e^/e^)/fiDc, where Dc = eBjm\yC is the cyclotron fre¬ 
quency, with mb being the electron band mass. For 
parameters appropriate for GaAs, namel y e ^ 12.5 
and mb ~ 0.067mg, we have A « 2.52 /y^ 5[T], which 
depends on the magnetic field and falls in the range 
A = 0.5 — 1.3 for the experimentally relevant fields 
5 = 4 — 25 T. For massless Dirac electrons in graphene, 
A = {e‘^fei)l{hvYfi) = ^jihevY) is essentially the 
graphene hue-structure constant, independent of mag¬ 
netic held. One obtains A ~ 2.2 for suspended graphene, 
A ~ 0.9 for graphene on Si 02 , A = 0.5 — 0.8 for graphene 
on Brvj™. 

We estimate the corrections due to LL mixing in the 
following manner. The ehect of LL mixing can be incor¬ 
porated into a LLL problem by modifying the interaction, 
by adding a “correction” term Vcorr that contains two-, 
three- and higher body interaction terms. The correc¬ 
tions to the hrst few psen dopoten tials of Korr have been 
evaluated in the literatur d^^ l ^^ l ^^ l in a perturbative treat¬ 
ment to linear order in A. We estimate the correction to 
the ground state energies by evaluating the expectation 
value of Vcorr with respect to the unperturbed ground 
state. We note that for the n = 0 graphene LL, no three 
body terms are induced. This method is expected to be 
valid for sufficiently small values of A. 

The technical details are as follows. (We include three- 
body interaction for completeness although it is not in¬ 
cluded in our calculations below.) For each FQHE state 


of interest, we evaluate the probability of occupation of 
various pair and triplet states, from which it is straight¬ 
forward to evaluate the correction to the energy using 
the two- and three-body pseudopotentials of Worr- To 
obtain the occupation amplitudes, we have computed by 
exact diagonalization a sequence of finite size Goulomb 
ground state vectors, labeled by N and 2Q = N/u — a 
(where a is the “shift” dependent on u and 7 ). For each 
vector, we then calculated the series of pair and triplet 
{K = 2 and 3) amplitudes for all possible 

pair and triplet spins (5 = 0 and lfori^ = 2;5 = l/2 
and 3/2 for if = 3) and the leading relative angular mo¬ 
menta (even m = 0, 2, ..., 8 for = 2 and 5 = 0; odd 
m = 1, 3, ..., 9 for = 2 and 5 = 1; m = 1, 2, 3 for 
K = 3 and s = 1/2; m = 3, 5, 6 , 7, 8, 9 for A" = 3 
and s = 3/2), as expectation values of the corresponding 
model K = 2 and 3 body pseudopotentials 
Owing to their regular size dependence, each amplitude 
was then reliably extrapolated by a linear regression as 
a function of 1/N to the limit of an infinite system to 

obtain Pi^\s,m) = limi/Ar->o-d, 7 ;V(®’ 

The maximum feasible dimension of about 4 x 10 ^ 
meant that for simple fractions we have data for many 
system sizes (e.g., for z/ = 2/3: N <23 for the polarized 
phase and A" < 14 for the unpolarized phase). However, 
for the more complex fractions such as u = 4/7 and 4/9 
we only have data for two sizes, with the smaller size 
suffering from the “aliasing” problem (e.g., for v = 4/9: 
A = 16 and 20 for the polarized phase, A = 10 and 
14 for the partially polarized phase, and A = 8 and 12 
for the unpolarized phase). The results for these systems 
are therefore less reliable. The fractions n/{2n ± 1) with 
n > 5 are not amenable to exact diagonalization studies. 

The extrapolated amplitudes were then convoluted 
with the effective pair and triplet LL mixing pseudopo¬ 
tentials Vc^r {s,m) derived in Ref. [69] (see Tables III 
and IV of that article) to give LL mixing corrections 
in the “linear” regime. In this regime, the ef¬ 
fect of LL mixing is estimated perturbatively, to the 
first order in A. Our estimates of LL mixing correc¬ 
tions to the ground state energy per particle are given per 
unit of A, separately for each filling factor and spin po- 
larizarion: Ae^,^/A = 

with the sums running over all spins and over the lead¬ 
ing angular momenta for which the pseudopotentials are 
availabl^^ (however, since m corresponds to an average 
K = 2 body distance or A = 3 body area, both K = 3 
sums (for 5 = 1/2 and 3/2) ought to be limited to the 
same rUmax = 3.). 

(It should also be mentioned that alternative meth¬ 
ods for including LL mixing are in principle possible. 
For example, one can attempt diagonalization of the A- 
electron Hamiltonian in an expanded Hilbert space, in¬ 
cluding cyclotron-excited con figura tions with some occu¬ 
pation of higher Landu levelsHowever, we have not 
found this method to be feasible for the present problem.) 

The modified critical values of /i: for A = 1 are shown 
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Filling factor u 

FIG. 3: (Color online) This figure compares the theoretical phase diagram with the experimentally measured one in graphene 
(taken from Feldman ei shown by green dots), GaAs (taken from Du ei shown by black squares) and AlAs systems 
(taken from Padmanabhan ei shown as red or blue triangles). The theoretical values (blue crosses) are from exact 

diagonalization at zero thickness and zero LL mixing. (We note that the blue crosses for the spin transitions at z/ = 4/7 and 
zy = 4/9 are obtained using results extrapolated from only two hnite systems.) For AlAs there is data available at two different 
densities: the blue downward triangle corresponds to density of 5.5 x 10^^ cm“^ while the red upward triangle corresponds to 
density of 5.0 x 10^^ cm“^. 


in the tables in the Appendix [A| for u = 2/3, 3/5, 4/7, 
2/5, 3/7 and 4/9. A comparison with the experimental 
values is shown in Fig. where in the experimental data 
was obtained by assuming ge = 6. A better agreement 
between the theoretical and experimental results can be 
obtained by choosing the value of ge = 16, but this value 
seems implausible. Therefore we come to the conclusion 
that theory substantially overestimates the effect of LL 
mixing. There may be several possible origins for this. 


First, the correction to the interaction has been calcu¬ 
lated to linear order in A perturbatively, and is thus valid 
only so long as the correction to the energy remains lin¬ 
ear in A. It is possible that A = 2.2 is outside the linear 
regime. We have also assumed that the wave functions 
themselves are not significantly modified by LL mixing. 
This should be the case for small A but may not be valid 
for A = 2.2. Finally, we have included corrections only 
for pseudopotentials up to a given relative angular mo¬ 
mentum. The quantitative errors due to such an ad hoc 
truncation are not known but may be significant. We 
believe that these comparisons bring out complications 
associated with the theoretical treatment of the quanti¬ 
tative effect of LL mixing. 


V. COMPOSITE FERMIONS CARRYING FOUR 
VORTICES 

It was predicted in Ref. [47] that the spin physics of the 
FQHE states at n/(4n+l), described in terms of compos¬ 
ite fermions carrying four vortices (^CFs), is qualitatively 
different from that at n/(2n ±1). Calculations based on 
the JK wave functions indicated that the fully spin polar¬ 
ized state (n, 0) at n/(4n + 1) is the ground state even at 
zero Zeeman energy, and consequently there are no spin- 
polarization phase transitions. The failure of the free-CF 
model was interpreted in Ref. |47| in terms of a Bloch fer¬ 
romagnetism for composite fermions, caused by a large 
exchange interaction that dominates their CF-cyclotron 
energy and favors the fully polarized state even in the 
absence of the Zeeman energy. 

We have seen above that the JK projection overesti¬ 
mates the energies of the non-fully spin polarized state 
by a larger amount than of the fully spin polarized states. 
One may therefore ask if the result in Ref. [43 is an arti¬ 
fact of the JK projection scheme. Furthermore, the spin 
physics at n/(4n — 1) has not been investigated so far. 
Are these states also always fully spin polarized? 

We have investigated these questions both by exact 
and CF diagonalizations for some of these fractions. For 
2/9 FQHE state, exact diagonalization results are incon- 
lusive but both the JK wave function and CE diagonal- 
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FIG. 4: (Color online) This figure compares the experimental phase diagram in graphene (taken from Feldman et aZ.^, shown 
by green dots) with the theoretical phase diagram (blue crosses) including corrections from LL mixing, assuming that the 
correction remains linear in A. The data of Feldman et al.^ assumes gt — 6.0. The theoretical estimates were obtained with 
A = 2.2, which corresponds to the suspended graphene samples of Feldman et al.^. (We note that the blue crosses for the spin 
transitions at z/ = 4/7 and z/ = 4/9 are obtained using results extrapolated from only two finite systems.) 


ization (see Fig. find that the fully spin polarized 
state has lower energy, thus confirming absence of any 
spin-polarization phase transition. The same is true for 
the 3/13 and 4/17 FQHE states wherein both the JK 
wave function and CF diagonalization find the state with 
larger spin polarization to be lower in energy than the 
state with a smaller spin polarization (see Fig. [^. We 
have also studied the first meaningful member of the Jain 
sequence n/(4n — 1), namely 2/7. Here, the calculations 
based on the JK wave functions as well as exact diago¬ 
nalization show that the spin singlet state has a slightly 
lower energy (see Fig. |^, giving /^ 2/7 = 0.0013(3) (ex¬ 
act). 

For the ^^CF states with p > 1, one can also perform 
the JK projection slightly differently as: 

(13) 

= j2P-2pLLL[$n]V2 (14) 


Because these wave functions apparently build better 
short distance correlations, one may expect them to 
have lower energies than their counterparts in Eqs. § 
and Contrary to this expectation we fin^that at 


u = nl{2pn + 1) the wave functions in Eq. 13 in gen¬ 
eral have slightly higher energies for finite systems than 
those in Eqs. [S, as seen in Eig. (We note that for 
states restricted to two A levels, the wave functions of 
Eq. and Eq. ^are identical.). In contrast, for states 
at z/ = n/{2pn — 1) we find that the ground state ener¬ 
gies obtained from the wave functions in Eq.f^are lower 


than those obtained from Eq. (see Eig. [^. Therefore, 
for states dXu = n/ (2pn+l) with p > 1 we quote energies 
obtained from Eqs. © while for states at = n/ {2pn — l) 
we quote energies obtained from Eq. 

VI. CONCLUSIONS 

The critical Zeeman energies where transitions be¬ 
tween differently spin polarized states occur are a direct 
measure of the energy difference between the states, and 
thus serve as a very sensitive test of the quantitative ac¬ 
curacy of the theory of the EQHE. These critical Zee- 
man energies are in general subject to corrections due 
to finite thickness of the quantum well and also Landau 
level mixing. In GaAs both effects are present, and it 
is not straightforward to disentangle their contributions, 
although progress has been made in experiments that 
study the effect systematically as a function of the quan¬ 
tum well widttP^. Because graphene has negligible finite 
thickness corrections, this gives an opportunity to obtain 
an accurate test of the CE theory, and also to gain insight 
into our understanding of the role of LL mixin^I^I^ 

We have evaluated accurate spin-polarization phase di¬ 
agram for the EQHE for an ideal two-dimensional system 
confined to the LLL with no LL mixing and no disorder. 
We have also evaluated corrections due to LL mixing, as¬ 
suming that these are linear in the parameter A. We find 
that the experimental results of Eeldman et are in 
excellent agreement with theory that neglects LL mix- 
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FIG. 5: (Color online) Thermodynamic extrapolation of the Coulomb ground state energies at v — 2/7 (top left panel), 2/9 
(top right panel), 3/13 (bottom left panel) and 4/17 (bottom right panel). The energies are obtained from exact diagonalization 
(“exact”), JK wave functions (“JK w.f.”) or CF diagonalization (“CFD”). 



FIG. 6: Comparison of the LLL Coulomb ground state en¬ 
ergies obtained from wave functions of Eq. (black squares) 
and Eq. (red squares) for fully polarized states at 3/13 
(left panel) and 4/17 (right panel). The former gives better 
energies and is used for the results in Tables [III] to IVIIll 



EIG. 7: Comparison of the LLL Coulomb ground state en¬ 
ergies obtained from wave functions of Eq. E and Eq. for 
states with reverse flux attachment. The latter gives better 
energies, and is used to obtain the energies shown in Tables 

[mlto lvlnl 


ing. Somewhat unexpectedly, if we include LL mixing in 
a linear approximation, the agreement becomes signifi¬ 
cantly worse, indicating that the amount of LL mixing 
in experiments is too large to be captured by a first or¬ 


der perturbative treatment. These results underscore our 
lack of a quantitative understanding of the effect of LL 
mixing on various quantities. 

We have shown that the critical Zeeman energies are 
well captured in terms of an effective mass model of com¬ 
posite fermions. We have shown that the CF states with 
Jain-Kamilla projection produce the correct energy or¬ 
dering for these states, and are fairly accurate for the 
“parallel flux attached” states at z/ = n/(2n + l), produc¬ 
ing the critical Zeeman energies with better than ^15% 
accuracy. In contrast, for the “reverse flux attached” 
states at = n/(2n — 1) the JK projection obtains the 
correct energy ordering of the states but obtains the crit¬ 
ical Zeeman energies that are off by a factor of two to 
three. For these states, the hard-core projection method 
is very accurate but difficult to implement for large sys¬ 
tems. 
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ly 

state 

N 

2Q 

S 

dimension 

2/3 

(2,0) 

28 

42 

14 

1,521,967,986 

2/3 

(1,1) 

14 

20 

0 

280,934,870 

2/5 

(2,0) 

18 

41 

9 

3,546,374,322 

2/5 

(1,1) 

12 

27 

0 

2,211,680,688 

3/5 

(3,0) 

24 

41 

12 

3,546,374,322 

3/5 

(2,1) 

14 

23 

3 

383,215,178 

3/7 

(3,0) 

18 

37 

9 

386,905,330 

3/7 

(2,1) 

11 

22 

2.5 

17,969,272 

4/7 

(4,0) 

20 

37 

10 

386,905,330 

4/7 

(3,1) 

14 

25 

5 

55,975,102 

4/7 

(2,2) 

12 

21 

0 

114,153,021 

4/9 

(4,0) 

20 

39 

10 

1,438,058,853 

4/9 

(3,1) 

14 

27 

5 

186,301,264 

4/9 

(2,2) 

12 

23 

0 

336,012,314 


TABLE II: Dimension of the Hilbert space of the largest sys¬ 
tems for which exact diagonalization results were obtained in 
this work. We show the dimension of states in the L = 0 
and relevant S sector for several values of (iV, 2Q) at various 
filling factors. 


Appendix A: Results 


In this Appendix, we give the results for the individ¬ 
ual systems used in obtaining the extrapolated energies 
shown in the main text. The tables below give ener¬ 
gies obtained from three methods. One is exact diago¬ 
nalization in the spherical geometry. These results are 
obtained by the Lanczos method. In Table |TI| we give 
the dimensions of the full lowest Landau level Hilbert 
space for the largest system size at various filling fac¬ 
tors considered in this work. These are the dimensions 
for the L = 0 and the relevant S sector. The results 
labeled “JK w.f’ are obtained from Jain’s CF wave func¬ 
tions using the JK projection method while those labeled 
“CFD” are obtained by the method of composite fermion 
diagonalizatiorP^. In Fig. |^we show the thermodynamic 
extrapolation of the ground state energies from finite size 
calculations of different spin polarized states at various 
filling factors in the lowest Landau level. Tables |III[ [IV| 
[V| |VI[ |VII| and VIII| show the thermodynamic energies 
obtained from these extrapolations and the critical Zee- 
man energies for the spin transitions. In these tables we 
also show the corrections obtained from LL mixing for 
the LL mixing parameter value A = 1. 


First one factorizes the Jastrow factor in Eq. (up to 
an overall sign), 




IT-'.. 




(Bl) 

(B2) 


Then each Ji is attached to the elements of a col¬ 
umn of the Slater determinant, and the projection 
is performed in each element individually. The S later 
determinant is composed of monopole harmonic^^lllll 
yQ,n,m{u^v)^ whcrc Q is the effective monopole strength 
for composite fermions, n is their A-level index and m 
is the value of the ^-component of the orbital angu¬ 
lar momentum operator. Projection of a single elec¬ 
tron wave function turns YQ^ri,m{ui,Vi) into an operator 
yQ,n,m{ui,Vi) that acts on the corresponding factor Ji. 

If 2p flux quanta are bound to each electron, p > 1, 
there one can follow two approaches. First, by Eqs. (13) 
and (14), 2p — 2 powers of J are moved outside of the 
LLL projection. Then YQ^n.m acts on Ji. Second, if by 
Eqs. ^ and ^ the complete Jastrow factor is within the 
scope of T^lll, Yq^n^m acts on Jf. For parallel flux at¬ 
tachment, the projected wave functions given in Refs.[45l 
|46] and [3] are applicable in both approaches. For reverse 
flux attachment {Q < 0), Davenport and SimorP^ gave 
an efficient method to obtain the projected wave func¬ 
tions. Because they implemented it explicitly only for 
p = 1, here we give, for completeness, LLL projection 
details for p = 2 for Q < 0. 

For p = 2 and Q < 0, we have 




,n,m 


n 

iUi,Vi)jf oc 

S = 0 



( 2|(5| + n \ 
\\Q\ TTl S J 




(B3) 


On the other hand. 


■^* = ( H ) E (B4) 

\j^ij t=o 


where 


Appendix B: Projections with reverse flux 
attachment 

We follow the JK method of projecting a CF wave 
function to the lowest Landau level. We use the 
spherical geometr}CHl Ui = cos{0i/2)e'^^Y‘^ and Vi = 
sm{ 0 i/ 2 )e~'^^Y‘^ are spinorial coordinates on the sphere. 


= 0,Ar-i 


Ui' 


Ui — l 



(B5) 


• • •, xm) 


0 otherwise 

(B6) 


Here, et,m are elementary symmetric polynomials. 
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v 

(2,0) 

(1,1) 



exact 

JK w.f. 

CFD 

exact 

JK w.f. 

CFD 

exact 

JK w.f. 

CFD 

2/3 

-0.51829(2) 

-0.5176(1) 


-0.52704(4) 

-0.5217(2) 

- 

0.0175(1) 

0.0065(17) 

- 

2/3 (with LL mixing) 

-0.61469(16) 

- 


-0.62969(10) 

- 

- 

0.0300(5) 

- 

- 

2/5 

-0.43298(3) 

-0.43277(2) 

-0.43287(2) 

-0.43935(1) 

-0.43839(2) 

-0.43902(3) 

0.0127(1) 

0.0113(1) 

0.0123(1) 

2/5 (with LL mixing) 

-0.46870(10) 

- 


-0.47932(17) 

- 

- 

0.0212(5) 

- 

- 

2/7 

-0.38185(9) 

-0.38140(6) 


-0.38249(8) 

-0.38188(6) 

- 

0.0013(3) 

0.0010(2) 

- 

2/9 

-0.34314(18) 

-0.34274(2) 

-0.34299(6) 

-0.34319(25) 

-0.34221(2) 

-0.34256(2) 

0.0001(9) 

- 

- 

2/11 

- 

-0.31331(2) 


- 

-0.31329(2) 


- 

-0.0002(3) 

- 

2/13 

- 

-0.29087(1) 

-0.29133(16) 

- 

-0.29019(1) 

-0.29047(3) 

- 


- 


TABLE III: Lowest Landau level Coulomb interaction energies (in units of jet) obtained from a thermodynamic extrapolation 
of results on the spherical geometry for various spin polarized states [denoted by (n^, n^)] states oX v — 2/ (4pdi 1). The energies 
are obtained from exact diagonalization (“exact”), JK wave functions (“JK w.L”) or CF diagonalization (“CFD”). Also shown 
are the critical Zeeman energies for spin transitions between the successive states. For Ezlie^ jet) > the state with larger 
polarization is favored over the one with smaller polarization. In this as well as the following tables, the energy corrections 
obtained from LL mixing are quoted for the value of the LL mixing parameter A=l. 


v 

(3,0) 

(2,1) 

K 


exact 

JK w.f. 

CFD 

exact 

JK w.f. 

CFD 

exact 

JK w.f. 

CFD 

3/5 

-0.49742(1) 

-0.4967(3) 


-0.50366(2) 

-0.4995(1) 

- 

0.0187(1) 

0.0081(12) 

- 

3/5 (with LL mixing) 

-0.57710(25) 

- 


-0.58873(46) 

- 

- 

0.0349(14) 


- 

3/7 

-0.44236(2) 

-0.4423(1) 

-0.44237(1) 

-0.44800(2) 

-0.44710(1) 

-0.44748(4) 

0.0167(2) 

0.0144(1) 

0.0154(2) 

3/7 (with LL mixing) 

-0.484663(13) 

- 


-0.49628(59) 

- 

- 

0.0349(14) 


- 

3/11 

- 

-0.3738(1) 


- 

-0.37352(7) 

- 

- 

- 

- 

3/13 

- 

-0.34839(3) 

-0.34844(1) 

- 

-0.34771(5) 

-0.34794(2) 

- 

- 

- 

3/17 

- 

-0.30924(7) 


- 

-0.30873(6) 

- 

- 

- 

- 


TABLE IV: Same as in Table III but for states at z/ = 3/(6p ± 1). 


v 

(4,0) 

(3,1) 

(2,2) 


exact 

JK w.f. 

CFD 

exact 

JK w.f. 

CFD 

exact 

JK w.f. 

CFD 

4/7 

-0.48842(0) 

-0.4875(7) 

- 

-0.49370(0) 

-0.4904(3) 

- 

-0.49495(0) 

-0.4908(2) 

- 

4/7 (with LL mixing) 

-0.56161(0) 

- 

- 

-0.57177(0) 


- 

-0.57432(0) 


- 

4/9 

-0.44771(0) 

-0.44750(1) 

-0.44770(10) 

-0.45241(0) 

-0.45155(1) 

-0.45184(3) 

-0.45382(0) 

-0.45275(2) 

-0.45288(6) 

4/9 (with LL mixing) 

-0.49303(0) 

- 

- 

-0.50256(0) 


- 

-0.50531(0) 


- 

4/17 

- 

-0.35123(1) 

-0.35125(2) 

- 

-0.35062(1) 

-0.35097(6) 

- 

-0.35053(1) 

-0.35051(1) 


TABLE V: Same as in Table III but for states at z/ = 4/(8p± 1). An underbar indicates that the thermodynamic extrapolation 
was done using only two systems, indicating that the results are less reliable. 


V 


1^2 


exact 

JK w.f. 

CFD 

exact 

JK w.f. 

CFD 

4/7 

0.0211(0) 

0.012(4) 

- 

0.0050(0) 

0.002(1) 

- 

4/7 (with LL mixing) 

0.0406(0) 

- 

- 

0.0102(0) 

- 

- 

4/9 

0.0188(0) 

0.0162(1) 

0.0166(5) 

0.0057(0) 

0.0048(1) 

0.0042(4) 

4/9 (with LL mixing) 

0.0381(0) 

- 

- 

0.0110(0) 

- 

- 


TABLE VI: The two critical Zeeman energies for spin transitions between the fully polarized and partially polarized states 
(z^i) and between the partially polarized and spin singlet states (z^ 2 ) at z/ = 4/(8p± 1). An underbar indicates that the 
thermodynamic extrapolation was done using only two systems, indicating that the results are less reliable. 
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V 

(5,0) 

(4,1) 

(3,2) 


f^2 


JK w.f. 

JK w.f. 

JK w.f. 

JK w.f. 

JK w.f. 

5/11 

-0.45080(1) 

-0.45429(2) 

-0.45563(4) 

0.0175(2) 

0.0067(3) 


TABLE VII: Same as in Table [nT| but for states at z/ = 5/(10p ± 1). Also shown are the two critical Zeeman energies for spin 
transitions between the fully polarized (5,0) and partially polarized (4,1) states (z^i) and between the partially polarized (4,1) 
and partially polarized (3,2) states (/^ 2 ). 


V 

(6,0) 

(5,1) 

(4,2) 

(3,3) 

Kx 

Z^2 

Z^3 


JK w.f. 

JK w.f. 

JK w.f. 

JK w.f. 

CED 

JK w.f. 

JK w.f. 

JK w.f. 

CED* 

6/13 

-0.45316(2) 

-0.45627(10) 

-0.45757(1) 

-0.45800(7) 

-0.45820(10) 

0.0186(7) 

0.0078(7) 

0.0026(5) 

0.0038(7) 


TABLE VIIL Same as in Table III but for states at z/ = 6/(12p ± 1). Also shown are the three critical Zeeman energies for 


spin transitions between the fully polarized and partially polarized (5,1) states (z^i); between the partially polarized (5,1) and 
partially polarized (4,2) states (z^ 2 ) and between the partially polarized (4,2) and spin-singlet states (z^a). The * for the CED 
value of k ,3 indicates that only the spin-singlet energy was calculated using CED; for the partially polarized (4,2) state, the 
zeroth order CE energy was used. 


S -0.540 

H 



• fully polarized (exact) 

• fully polarized (JK w.f) 

• spin singlet (exact) 

• spin singlet (JK w.f) 

i/ = 2/3 



0.00 0.05 


0.15 

1/N 



-0.500 


• fully polarized (exact) 

• fully polarized (JK w.f) 

• partially polarized (exact) 

• partially polarized (JK w.f) 

-0.430 

• fully polarized (exact) 

• fully polarized (JK w.f) 


-0.480 

- u = 4/7 ^-- 




fully polarized (CFD) 

• partially polarized (exact) 

• partially polarized (JK w.f) 

• partially polarized (CFD) 






5-0.435 


^-0.490 


/:/-0.520 




A 


0) 


bC 

^ -0.440 

& 


^-0.500 

c 

H 

• (4,0) (exact) 

• (4,0) (JK w.f) 

• (3,1) (exact) 


-0.540 

z/ = 3/5 ^ 

-0.445 


I' = 3/7 

-0.510 

• (3,1) (JK w.f) 

• (2,2) (exact) 







(2,2) (JK w.f) 



0.10 0.15 

1/N 


0.10 

1/N 


1/N 



EIG. 8: Extrapolation of the lowest Landau level Coulomb ground state energy to the thermodynamic limit for different spin 
polarized states at various filling factors. 
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Thus, when evaluating Eq. |B3[ the key step is 


where Pr^N is the power-sum polynomial defined as: 


/N-1 \ /N-1 

E(-l) 




, t=0 


^ r=0 


( \ N-1 N-1 

j^i J t=0 r=0 

= fn“d 

\j7^i J t,r=0 
{2{N -1) - t- 

{2{N — 1) — t — r — {\Q\ — m -\- n — s))! 


{t + r - (IQI + m + s))!' 


(B7) 


N 

Pr,N{xi,X2,-■ ■ ,Xn) = '^xl (Bll) 

i=l 

We also note another iterative identity of the symmetric 
polynomials: 


em,N-l{xi,X2, • • • ,Xj^i • • • ,Xn) = em,N{xi,X2, * * * ,Xn) 
1,V—1 (^1: ^2 5 ' ' ' ^ Xjj^i ' ' ' , Xw ) 

The above two identities can be used in conjunction with 
each other to create an efficient routine to store the com¬ 
plete set of symmetric polynomials ej. A word of caution 
is due here: the above quantities tend to suffer from nu¬ 
merical precision errors. To get around this problem, we 
store all the numerical values to high precision. 

We also give the projection formula for p = 1 correcting 
a typo in Ref. |62j 


The summation in Eq. |B7| must be restricted as 

t + r > |Q| + m + 5 , (B 8 ) 

t + r < 2(At — 1) — (IQI — m + n — s). (B9) 

The elementary symmetric polynomials stated above 
can be calculated iteratively using the following Newton’s 
identit}^^ 

^ m 

em,N{xi,X2,-■ ■ ,Xn) = —N{-'i-Y''~^Pr,N{xi,X2,---,XN) 

m 


n 

^Q,n,rnij^i-)Xi)Ji OC ^ ^ 1) 


j^i s =0 
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E (-1) 

t — I Q I 


tA 


t\ 


' [t - (IQI + m + s)]! 


n\ / 2\Q\+n 
s) \|Q| +m + s 


t-(|Q|+m)^ 


[AT - 1 -t]! 


N—l—t—(\Q\—m) 


[N — 1 — t — (IQI — m ^ n — 5)]! 


(B12) 


^ r,V(:^l 5 ^ 2 : ■ ' ■ 5 ^V) (BIO) 
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